Evolutionary biology of Gondwanan non-biting midges (Diptera: Chironomidae)

The potential restriction to effective dispersal and gene flow caused by habitat fragmentation can apply to multiple levels of evolutionary scale; from the fragmentation of ancient supercontinents driving diversification and speciation on disjunct landmasses, to the isolation of proximate populations as a result of their inability to cross intervening unsuitable habitat. Investigating the role of habitat fragmentation in driving diversity within and among taxa can thus include inferences of phylogenetic relationships among taxa, assessments of intraspecific phylogeographic structure and analyses of gene flow among neighbouring populations. The proposed Gondwanan clade within the chironomid (non-biting midge) subfamily Orthocladiinae (Diptera: Chironomidae) represents a model system for investigating the role that population fragmentation and isolation has played at different evolutionary scales. A pilot study by Krosch et al (2009) indentified several highly divergent lineages restricted to ancient rainforest refugia and limited gene flow among proximate sites within a refuge for one member of this clade, Echinocladius martini Cranston. This study provided a framework for investigating the evolutionary history of this taxon and its relatives more thoroughly. Populations of E. martini were sampled in the Paluma bioregion of northeast Queensland to investigate patterns of fine-scale within- and among-stream dispersal and gene flow within a refuge more rigorously. Data was incorporated from Krosch et al (2009) and additional sites were sampled up- and downstream of the original sites. Analyses of genetic structure revealed strong natal site fidelity and high genetic structure among geographically proximate streams. Little evidence was found for regular headwater exchange among upstream sites, but there was distinct evidence for rare adult flight among sites on separate stream reaches. Overall, however, the distribution of shared haplotypes implied that both larval and adult dispersal was largely limited to the natal stream channel. Patterns of regional phylogeographic structure were examined in two related austral orthoclad taxa – Naonella forsythi Boothroyd from New Zealand and Ferringtonia patagonica Sæther and Andersen from southern South America – to provide a comparison with patterns revealed in their close relative E. martini. Both taxa inhabit tectonically active areas of the southern hemisphere that have also experienced several glaciation events throughout the Plio-Pleistocene that are thought to have affected population structure dramatically in many taxa. Four highly divergent lineages estimated to have diverged since the late Miocene were revealed in each taxon, mirroring patterns in E. martini; however, there was no evidence for local geographical endemism, implying substantial range expansion post-diversification. The differences in pattern evident among the three related taxa were suggested to have been influenced by variation in the responses of closed forest habitat to climatic fluctuations during interglacial periods across the three landmasses. Phylogeographic structure in E. martini was resolved at a continental scale by expanding upon the sampling design of Krosch et al (2009) to encompass populations in southeast Queensland, New South Wales and Victoria. Patterns of phylogeographic structure were consistent with expectations and several previously unrecognised lineages were revealed from central- and southern Australia that were geographically endemic to closed forest refugia. Estimated divergence times were congruent with the timing of Plio-Pleistocene rainforest contractions across the east coast of Australia. This suggested that dispersal and gene flow of E. martini among isolated refugia was highly restricted and that this taxon was susceptible to the impacts of habitat change. Broader phylogenetic relationships among taxa considered to be members of this Gondwanan orthoclad group were resolved in order to test expected patterns of evolutionary affinities across the austral continents. The inferred phylogeny and estimated divergence times did not accord with expected patterns based on the geological sequence of break-up of the Gondwanan supercontinent and implied instead several transoceanic dispersal events post-vicariance. Difficulties in appropriate taxonomic sampling and accurate calibration of molecular phylogenies notwithstanding, the sampling regime implemented in the current study has been the most intensive yet performed for austral members of the Orthocladiinae and unsurprisingly has revealed both novel taxa and phylogenetic relationships within and among described genera. Several novel associations between life stages are made here for both described and previously unknown taxa. Investigating evolutionary relationships within and among members of this clade of proposed Gondwanan orthoclad taxa has demonstrated that a complex interaction between historical population fragmentation and dispersal at several levels of evolutionary scale has been important in driving diversification in this group. While interruptions to migration, colonisation and gene flow driven by population fragmentation have clearly contributed to the development and maintenance of much of the diversity present in this group, long-distance dispersal has also played a role in influencing diversification of continental biotas and facilitating gene flow among disjunct populations.

Dispersal, sex ad clonality in the marine environment: population genetic structure of the seagrass Cymodocea nodosa on Mediterranean and Atlantic coasts

Dissertação de dout. em Ecologia, Faculdade de Ciências do Mar e do Ambiente, Univ. do Algarve, 2005

Comparison of intraspecific genetic structure among related chironomids (Diptera) from New Zealand and Patagonia : disparity between potential and realized dispersal

Population genetic studies of freshwater invertebrate taxa in New Zealand and South America are currently few despite the geologically and climatically dynamic histories of these regions. The focus of our study was a comparison of the influence on realized dispersal of 2 closely related nonbiting midges (Chironomidae) of population fragmentation on these separated austral land masses. We used a 734-base pair (bp) fragment of cytochrome c oxidase subunit I (COI) to investigate intraspecific genetic structure in Naonella forsythi Boothroyd in New Zealand and Ferringtonia patagonica Edwards in Patagonia. We proposed hypotheses about their potential dispersal and, hence, expected patterns of genetic structure in these 2 species based on published patterns for the closely related Australian taxon Echinocladius martini Cranston. Genetic structure revealed for both N. forsythi and F. patagonica was characterized by several highly divergent (2.0–10.5%) lineages of late Miocene–Pliocene age within each taxon that were not geographically localized. Many were distributed widely. This pattern differed greatly from population structure in E. martini, which was typified by much greater endemicity of divergent genetic lineages. Nevertheless, diversification of lineages in all 3 taxa appeared to be temporally congruent with the onset of late Miocene glaciations in the southern hemisphere that may have driven fragmentation of suitable habitat, promoting isolation of populations and divergence in allopatry. We argue that differences in realized dispersal post-isolation may be the result of differing availability of suitable habitat in interglacial periods.

Population genetic structure and patterns of dispersal in the Giant Long-Armed Prawn, Macrobrachium lar (Fabricius, 1798) (Decapoda : Palaemonidae)

The Giant Long-Armed Prawn, Macrobrachium lar is a freshwater species native to the Indo-Pacific. M. lar has a long-lived, passive, pelagic marine larval stage where larvae need to colonise freshwater within three months to complete their development. Dispersal is likely to be influenced by the extensive distances larvae must transit between small oceanic islands to find suitable freshwater habitat, and by prevailing east to west wind and ocean currents in the southern Pacific Ocean. Thus, both intrinsic and extrinsic factors are likely to influence wild population structure in this species. The present study sought to define the contemporary broad and fine-scale population genetic structure of Macrobrachium lar in the south-western Pacific Ocean. Three polymorphic microsatellite loci were used to assess patterns of genetic variation within and among 19 wild adult sample sites. Statistical procedures that partition variation implied that at both spatial scales, essentially all variation was present within sample sites and differentiation among sites was low. Any differentiation observed also was not correlated with geographical distance. Statistical approaches that measure genetic distance, at the broad-scale, showed that all south-western Pacific Islands were essentially homogeneous, with the exception of a well supported divergent Cook Islands group. These findings are likely the result of some combination of factors that may include the potential for allelic homoplasy, through to the effects of sampling regime. Based on the findings, there is most likely a divergent M. lar Cook Islands clade in the south-western Pacific Ocean, resulting from prevailing ocean currents. Confirmation of this pattern will require a more detailed analysis of nDNA variation using a larger number of loci and, where possible, use of larger population sizes.

Population genetic structure and dispersal in white-lipped peccaries (Tayassu pecari) from the Brazilian Pantanal

In many mammals social organization promotes genetic structuring, which can be influenced by the dispersal pattern of the species. We analyzed the population genetic structure and dispersal of white-lipped peccaries (Tayassu pecan) from the Pantanal, Brazil. We genotyped 100 individuals at 7 microsatellite loci from 2 adjacent locations with no obvious geographic barrier between them. We found a significant but low F(ST) value, and the Bayesian analysis indicated a unique cluster. No significant differences were observed between mean assignment indices of resident males and females from both locations, and the probability of being born at the location sampled of > 30% of the individuals analyzed was lower than average. Mean relatedness between resident female, male, and opposite-sex pairs was not statistically different in both locations. These results suggest a low degree of genetic differentiation between the locations analyzed, and dispersal by both sexes (contrary to the predicted male-biased dispersal of most mammalian species).

Population genetic structure and dispersal in white-lipped peccaries (Tayassu pecari) from the Brazilian Pantanal

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Genetic effects of chronic habitat fragmentation on tree species: the case of Sorbus aucuparia in a deforested Scottish landscape

Sustainable forest restoration and management practices require a thorough understanding of the influence that habitat fragmentation has on the processes shaping genetic variation and its distribution in tree populations. We quantified genetic variation at isozyme markers and chloroplast DNA (cpDNA), analysed by polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) in severely fragmented populations of Sorbus aucuparia (Rosaceae) in a single catchment (Moffat) in southern Scotland. Remnants maintain surprisingly high levels of gene diversity (H-E) for isozymes (H-E = 0.195) and cpDNA markers (H-E = 0.490). Estimates are very similar to those from non-fragmented populations in continental Europe, even though the latter were sampled over a much larger spatial scale. Overall, no genetic bottleneck or departures from random mating were detected in the Moffat fragments. However, genetic differentiation among remnants was detected for both types of marker (isozymes Theta(n) = 0.043, cpDNA Theta(c) = 0.131; G-test, P-value < 0.001). In this self-incompatible, insect-pollinated, bird-dispersed tree species, the estimated ratio of pollen flow to seed flow between fragments is close to 1 (r = 1.36). Reduced pollen-mediated gene flow is a likely consequence of habitat fragmentation, but effective seed dispersal by birds is probably helping to maintain high levels of genetic diversity within remnants and reduce genetic differentiation between them.

Genetic dynamics across early life stages in the tropical tree Prunus africana (Rosaceae)

In many plant species, the genetic template of early life-stages is formed by animal-mediated pollination and seed dispersal and has profound impact on further recruitment and population dynamics. Understanding the impact of pollination and seed dispersal on genetic patterns is a central issue in plant population biology. In my thesis, I investigated (i) contemporary dispersal and gene flow distances as well as (ii) genetic diversity and spatial genetic structure (SGS) across subsequent recruitment stages in a population of the animal-pollinated and dispersed tree Prunus africana in Kakamega Forest, West Kenya. Using microsatellite markers and parentage analyses, I inferred distances of pollen dispersal (father-to-mother), seed dispersal/maternal gene flow (mother-to-offspring) as well as paternal gene flow (father-to-offspring) for four early life stages of the species (seeds and fruits, current year seedlings, seedlings ≤ 3yr, seedlings > 3yr). Distances of pollen and seed dispersal as well as paternal gene flow were significantly shorter than expected from the spatial arrangement of trees and sampling plots. They were not affected by the density of conspecific trees in the surrounding. At the propagule stage, mean pollen dispersal distances were considerably (23-fold) longer than seed dispersal distances, and paternal gene flow distances exceeded maternal gene flow by a factor of 25. Seed dispersal distances were remarkably restricted, potentially leading to a strong initial SGS. The initial genetic template created by pollination and seed dispersal was extensively altered during later recruitment stages. Potential Janzen-Connell effects led to markedly increasing distances between offspring and both parental trees in older life stages. This showed that distance and density-dependent mortality factors are not exclusively related to the mother tree, but also to the father. Across subsequent recruitment stages, the pollen to seed dispersal ratio and the paternal to maternal gene flow ratio dropped to 2.1 and 3.4, respectively, in seedlings > 3yr. The relative changes in effective pollen dispersal, seed dispersal, and paternal gene flow distances across recruitment stages elucidate the mechanisms affecting the contribution of the two processes pollen and seed dispersal to overall gene flow. Using the same six microsatellite loci, I analyzed genetic diversity and SGS across five life stages, from seed rain to adults. Levels of genetic diversity within the studied P. africana population were comparable to other Prunus species and did not vary across life stages. In congruence with the short seed dispersal distances, I found significant SGS in all life stages. SGS decreased from seed and early seedling stages to older juvenile stages, and it was higher in adults than in late juveniles of the next generation. A comparison of the data with direct assessments of contemporary gene flow patterns indicate that distance- or density-dependent mortality, potentially due to Janzen-Connell effects, led to the initial decrease in SGS. Intergeneration variation in SGS could have been driven by variation in demographic processes, the effect of overlapping generations, and local selection processes. Overall, my study showed that complex sequential processes during recruitment contribute to the spatial genetic structure of tree populations. It highlights the importance of a multistage perspective for a comprehensive understanding of the impact of animal-mediated pollen and seed dispersal on spatial population dynamics and genetic patterns of trees.

The development of 10 novel polymorphic microsatellite markers through next generation sequencing and a preliminary population genetic analysis for the endangered Glenelg spiny crayfish, Euastacus bispinosus

The Glenelg spiny crayfish, Euastacus bispinosus, is an iconic freshwater invertebrate of south eastern Australia and listed as 'endangered' under the Environment Protection and Biodiversity Conservation Act 1999, and 'vulnerable' under the International Union for Conservation of Nature's Red List. The species has suffered major population declines as a result of over-fishing, low environmental flows, the introduction of invasive fish species and habitat degradation. In order to develop an effective conservation strategy, patterns of gene flow, genetic structure and genetic diversity across the species distribution need to be clearly understood. In this study we develop a suite of polymorphic microsatellite markers by next generation sequencing. A total of 15 polymorphic loci were identified and 10 characterized using 22 individuals from the lower Glenelg River. We observed low to moderate genetic variation across most loci (mean number of alleles per locus = 2.80; mean expected heterozygosity = 0.36) with no evidence of individual loci deviating significantly from Hardy-Weinberg equilibrium. Marker independence was confirmed with tests for linkage disequilibrium, and analyses indicated no evidence of null alleles across loci. Individuals from two additional sites (Crawford River, Victoria; Ewens Ponds Conservation Park, South Australia) were genotyped at all 10 loci and a preliminary investigation of genetic diversity and population structure was undertaken. Analyses indicate high levels of genetic differentiation among sample locations (F ST = 0.49), while the Ewens Ponds population is genetically homogeneous, indicating a likely small founder group and ongoing inbreeding. Management actions will be needed to restore genetic diversity in this and possibly other at risk populations. These markers will provide a valuable resource for future population genetic assessments so that an effective framework can be developed for implementing conservation strategies for E. bispinosus.

Spatial population genetic structure reveals strong natal site fidelity in Echinocladius martini (Diptera: Chironomidae) in northeast Queensland, Australia

1. A diverse array of patterns has been reported regarding the spatial extent of population genetic structure and effective dispersal in freshwater macroinvertebrates. In river systems, the movements of many taxa can be restricted to varying degrees by the natural stream channel hierarchy. 2. In this study, we sampled populations of the non-biting freshwater midge Echinocladius martini in the Paluma bioregion of tropical northeast Queensland to investigate fine scale patterns of within- and among-stream dispersal and gene flow within a purported historical refuge. We amplified a 639 bp fragment of mitochondrial COI and analysed genetic structure using pairwise ΦST, hierarchical AMOVA, Mantel tests and a parsimony network. Genetic variation was partitioned among stream sections using Streamtree to investigate the effect of potential instream dispersal barriers. 3. The data revealed strong natal site fidelity and significant differentiation among neighbouring, geographically proximate streams. We found evidence for only episodic adult flight among sites on separate stream reaches. Overall, however, our data suggested that both larval and adult dispersal was largely limited to within a stream channel. 4. This may arise from a combination of the high density of riparian vegetation physically restricting dispersal and from the joint effects of habitat stability and large population sizes. Together these may mitigate the requirement for movement among streams to avoid inbreeding and local extinction due to habitat change and may thus enable persistence of upstream populations in the absence of regular compensatory upstream flight. Taken together, these data suggest that dispersal of E. martini is highly restricted, to the scale of only a few kilometres, and hence occurs predominantly within the natal stream.

Phylogeography of Echinocladius martini Cranston (Diptera: Chironomidae) in closed forest streams of eastern Australia

The eastern Australian rainforests have experienced several cycles of range contraction and expansion since the late Miocene that are closely correlated with global glaciation events. Together with ongoing aridification of the continent, this has resulted in current distributions of native closed forest that are highly fragmented along the east coast. Several closed forest endemic taxa exhibit patterns of population genetic structure that are congruent with historical isolation of populations in discrete refugia and reflect evolutionary histories dramatically affected by vicariance. Currently, limited data are available regarding the impact of these past climatic fluctuations on freshwater invertebrate taxa. The non-biting midge species Echinocladius martini Cranston is distributed along the east coast and inhabits predominantly montane streams in closed forest habitat. Phylogeographic structure in E. martini was resolved here at a continental scale by incorporating data from a previous pilot study and expanding the sampling design to encompass populations in the Wet Tropics of north-eastern Queensland, south-east Queensland, New South Wales and Victoria. Patterns of phylogeographic structure revealed several deeply divergent mitochondrial lineages from central and south-eastern Australia that were previously unrecognised and were geographically endemic to closed forest refugia. Estimated divergence times were congruent with late Miocene onset of rainforest contractions across the east coast of Australia. This suggested that dispersal and gene flow among E. martini populations isolated in refugia has been highly restricted historically. Moreover, these data imply, in contrast to existing preconceptions about freshwater invertebrates, that this taxon may be acutely susceptible to habitat fragmentation.

Systematics and biogeography of the Gondwanan Orthocladiinae (Diptera: Chironomidae)

Restrictions to effective dispersal and gene flow caused by the fragmentation of ancient supercontinents are considered to have driven diversification and speciation on disjunct landmasses globally. Investigating the role that these processes have played in the development of diversity within and among taxa is crucial to understanding the origins and evolution of regional biotas. Within the chironomid (non-biting midge) subfamily Orthocladiinae (Diptera: Chironomidae), a group of genera that are distributed across the austral continents (Australia, New Zealand, South America) have been proposed to represent a relict Gondwanan clade. We used a molecular approach to resolve relationships among taxa with the aim to determine the relative roles that vicariance and dispersal may have played in the evolution of this group. Continental biotas did not form monophyletic groups, in accordance with expectations given existing morphological evidence. Patterns of phylogenetic relationships among taxa did not accord with expected patterns based on the geological sequence of break-up of the Gondwanan supercontinent. Likewise, divergence time estimates, particularly for New Zealand taxa, largely post-dated continental fragmentation and implied instead that several transoceanic dispersal events may have occurred post-vicariance. Passive dispersal of gravid female chironomid adults is the most likely mechanism for transoceanic movement, potentially facilitated by West Wind Drift or anti-cyclone fronts. Estimated timings of divergence among Australian and South American Botryocladius, on the other hand, were congruent with the proposed ages of separation of the two continents from Antarctica. Taken together, these data suggest that a complex relationship between both vicariance and dispersal may explain the evolution of this group. The sampling regime we implemented here was the most intensive yet performed for austral members of the Orthocladiinae and unsurprisingly revealed several novel taxa that will require formal description.

Not drowning, (hand)waving? Molecular phylogenetics, biogeography and evolutionary tempo of the ‘Gondwanan’ midge Stictocladius Edwards (Diptera: Chironomidae)

Many insect clades, especially within the Diptera (true flies), have been considered classically ‘Gondwanan’, with an inference that distributions derive from vicariance of the southern continents. Assessing the role that vicariance has played in the evolution of austral taxa requires testing the location and tempo of diversification and speciation against the well-established predictions of fragmentation of the ancient super-continent. Several early (anecdotal) hypotheses that current austral distributions originate from the breakup of Gondwana derive from studies of taxa within the family Chironomidae (non-biting midges). With the advent of molecular phylogenetics and biogeographic analytical software, these studies have been revisited and expanded to test such conclusions better. Here we studied the midge genus Stictocladius Edwards, from the subfamily Orthocladiinae, which contains austral-distributed clades that match vicariance-based expectations. We resolve several issues of systematic relationships among morphological species and reveal cryptic diversity within many taxa. Time-calibrated phylogenetic relationships among taxa accorded partially with the predicted tempo from geology. For these apparently vagile insects, vicariance-dated patterns persist for South America and Australia. However, as often found, divergence time estimates for New Zealand at c. 50 mya post-date separation of Zealandia from Antarctica and the remainder of Gondwana, but predate the proposed Oligocene ‘drowning’ of these islands. We detail other such ‘anomalous’ dates and suggest a single common explanation rather than stochastic processes. This could involve synchronous establishment following recovery from ‘drowning’ and/or deleteriously warming associated with the mid-Eocene climatic optimum (hence ‘waving’, which refers to cycles of drowning events) plus new availability of topography providing of cool running waters, or all these factors in combination. Alternatively a vicariance explanation remains available, given the uncertain duration of connectivity of Zealandia to Australia–Antarctic–South America via the Lord Howe and Norfolk ridges into the Eocene.

Drift dispersal in central Southern Ontario /

Two tills are readily identi-f i able in central Southern Ontario, a very stony, loose deposit o-f variable matrix (Dummer till) and a moderately stony, fissile and compact deposit that is more homogeneous (drumlinized till). The quantity o-f Precambr i an, Paleozoic and Shadow Lake Formation (Paleozoic) rock types were determined and corresponding isopleth maps drawn. The changes in lithology content occurred in the direction o-f transport, there-fore, compositional isopleths o-f till may be considered equipotential lines for the reconstruction of glacier flow paths. Areal gradations of drift lithology indicated that the prime agents of dispersal were ice and glacial meltwaters. The down-ice abundance trend of till components indicated a dispersal pattern showing the concentration of a given lithology type peaking within a few kilometres of the source followed by a rapid decline and thereafter, a more gradual decrease with increasing distance. Within the esker deposits, igneous rocks may form the major component and can extend further onto the limestone plain than in the adjacent till. Evidence is presented that indicates the "style" of dispersal was one in which glacial ice may have been strongly influenced by local bedrock topography and the regional structural trends. The ice tended to follow pre-existing valleys and lows, depositing till composed mainly of local bedrock. Gradations in Paleozoic clast content showed that the local bedrock lithology became the primary till component within 3 km of down-ice transport. Evidence is presented that indicated the last glaciation may have occurred as a relatively thin ice mass, followed by stagnation and recession. No evidence of a lateglacial re-advance was found within the study area. Because of the lack of a contact between the Dummer and drumlinized till, and because of results showing gradation of the Dummer till into the drumlinized till (as indicated by lithology content and grain size), it is suggested that no re-advance occurred.

Population genetic structure, gene flow and sex-biased dispersal in frillneck lizards (Chlamydosaurus kingii)

By using both mitochondrial and nuclear multiloci markers, we explored population genetic structure, gene flow and sex-specific dispersal of frillneck lizards (Chlamydosaurus kingii) sampled at three locations, separated by 10 to 50 km, in a homogenous savannah woodland in tropical Australia. Apart from a recombinant lizard, the mitochondrial analyses revealed two nonoverlapping haplotypes/populations, while the nuclear markers showed that the frillneck lizards represented three separate clusters/populations. Due to the small population size of the mtDNA, fixation may occur via founder effects and/or drift. We therefore suggest that either of these two processes, or a combination of the two, are the most likely causes of the discordant results obtained from the mitochondrial and the nuclear markers. In contrast to the nonoverlapping mitochondrial haplotypes, in 12 out of 74 lizards, mixed nuclear genotypes were observed, hence revealing a limited nuclear gene flow. Although gene flow should ultimately result in a blending of the populations, we propose that the distinct nuclear population structure is maintained by frequent fires resulting in local bottlenecks, and concomitant spatial separation of the frillneck lizard populations. Limited mark-recapture data and the difference in distribution of the mitochondrial and nuclear markers suggest that the mixed nuclear genotypes were caused by juvenile male-biased dispersal.